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Abstract 

We present in this continuation paper a new axiomatic derivation 
of the Schrodinger equation from three basic postulates. This new 
derivation sheds some light on the thermodynamic character of the 
quantum formalism. We also show the formal connection between 
this derivation and the one previously done by other means. Some 
considerations about metaestability are also drawn. We return to an 
example previously developed to show how the connection between 
both derivations works. 
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1 Introduction 



In the first paper of this series we showed how to mathematically derive 
the Schrodinger equation from three basic postulates. At that time we have 
made use of an infinitesimal Fourier transformation called by ourselves the 
Wigner-Moyal infinitesimal transformation. 

The use of this transformation, although very practical for formal devel- 
opments H2fl-|14j|) implies in some hiding of the underlying physics. 



In this continuation paper we will show how to derive the Schrodinger 
equation from more physically sounded axioms. This will not only fix the 
derivation itself, but also improve our understanding of the physics involved, 
which is our primary goal. 

This will be done in the second section. We will show that from three 
clear physical postulates it is possible to derive the Schrodinger equation. 

Then, in section three, we will present the formal (and physical) connec- 
tions between the present derivation and the one made in the first paper of 
this series [Q. 

Section four will be devoted to a brief consideration about the dynamics 
of the ensemble statistics. We will then show that there is a corresponding 
Hamilton equation for the ensemble dynamics. 

In section five, we show how our present formulation of the quantum me- 
chanical postulates enables us to treat metastablility behavior on qualitative 
grounds, by means of the entropy concept. 

The sixth section will then use the Infinitesimal Transformation formal- 
ism already developed [|I|] to deal with the metastability problem in a totally 
formal fashion. 

The last section is devoted to our conclusions. 

In the appendix, the basic concepts about the present Schrodinger equa- 
tion derivation, together with its connection with the previous derivation, 
will be exemplified with the concrete problem of a Boltzmann distribution. 
We will then show that this probability density has the related probability 
amplitudes (wave functions) that are solutions of a Schrondiger equation. A 
special application of this last result is then shown for the harmonic oscilla- 
tor. This example has already been developed in a previous paper ]TJ| but 
now will be used to show that it is a particular case of a much wider physical 
framework. 
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2 Derivation 



We begin with the Liouville equation 

dF(x,t;t) p dF(x,p;t) dV(x) dF(x,p;t) 



0. 



1) 



dt m dx dx dp 

Integrating this equation with respect to p and using the identities 

J F(x,p;t)dp = p(x;t) ; J pF(x,p;t)dp = p(x; t) p(x; t) , (2) 

where p(x; t) is the probability density on the configuration space and p(x; t) 

(3) 



is the so called macroscopic momentum [15|, we finally get 



representing a continuity equation for the probability density p(x;t). 

Now, we may multiply the Liouville equation (|I|) by p and integrate with 
respect to this variable to get, using 



p 2 F(x,p;t)dp = M 2 (x;t), 



(4) 



the equation 



d . , w X1 1 dM 2 (x;t) dV(x) . s 

s W; t)p(r, t)] + --^i + -gUpfc t) = 0. (5) 

Using the equation @ into equation (§) we get, after some straightfor- 
ward calculations, the equation 



m 9x 



M 2 (x; t) — p 2 (x; t)p(x; t) 



+ 



+ p{x;t) 



dp(x;t) d fp 2 {x; t)\ dV(x) 
dt dx \ 2m I dx 



(6) 



The first term to the left may be written as 



M 2 (x;t) — p 2 (x; t) p(x; t) = / p 2 —p 2 (x;t) F(x,p;t)dp 
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[p — p(x; t)] 2 F(x, p; t)dp. 



We now put 



(5p) 2 p(x;t) = / \p — p(x;t)\ F(x,p;t)dp 



(7) 



(8) 



and try to find a functional expression for (5p) 2 . 

To achieve this goal, we begin considering the entropy of this isolated 
system as given by 

S(x;t) = fcln[n(x;t)], (9) 

where Q(x; t) represents the system accessible states when the position x 
varies between x and x + Sx. If this position is free to vary ||16|| , then the 
equal a priori probability postulate says that 



p{x;t) (xQ{x;t) = e s ^ t)/k . 



(10) 



Now, we may ask which probability density p eq (x;t), defined upon the 
configuration space, is related with the thermodynamic equilibrium of the 
system. We are thus interested in the functional variation Ap(x; t) repre- 
senting densities variations induced by the fluctuations [|HJ . 

Then we expand the entropy around a thermodynamic equilibrium con- 
figuration, where the density is given by p eq (x;t), as 



S(x ± 5x; t) = S eq (x; t) + 



1 f d 2 S eq (x;tY 

2 \ dx 2 



{5xf + .. , 



(11) 



5x=0 



where S eq (x; t) stands for the thermodynamic equilibrium configuration en- 
tropy, given by 5x = 0, and where we have already used the fact that, for 
this configuration the entropy shall be a maximum and so 



f dS eq (x;t) ' 
\ dx 

We thus use the fact that 

p eq (x;t) = • 

to get 

p(x,5x;t) = p eq (x;t)exp - 



0. 



<5x=0 



D S eq (x;t)/k 



1 

2k 



d 2 S eq (x;t) 



dx 2 



(Sx) 



(12) 

(13) 
(14) 
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where we have already used the fact that the second order derivative of 
the entropy must be negative near an equilibrium point, together with the 
equation (|T2{). 

Thus, at each point x, the probability distribution, with respect to the 
small displacements Sx(x; t), is a gaussian and is related with the probability 
of having a fluctuation Ap, on the thermodynamic equilibrium density, owing 
to a fluctuation 5x. Equation (|14]) then guarantees that the system will tend 
to return to its equilibrium distribution represented by the probability density 

Peq t) ■ 

In this case, using the equation (]T4|), the mean quadratic displacements 
related with the fluctuations are given by 



(5x) 2 

where we have made 



f^(5x) 2 e-^ 2 d(5x) 
!+™e-^) 2 d{8x) 



1 

2? 



f 

27 



d 2 S eq (x;t) 



dx 2 



(15) 



(16) 



We have, a priori, no relation between these displacement fluctuations 
and those related with the momenta. We now impose the restriction that, in 
this thermodynamic equilibrium situation, we must have 



(dp) 2 (5x) 



4 ' 



We then get 



(5p) 2 = -— 



d 2 ln[p(x;t)] 



dx 2 



With this last result we may write 



n 2 



(Sp) 2 p(x;t) = -—p{x;t) 



d 2 ln[p(x; t)] 

dx 2 



Substituting this expression into equation @ and putting 

ds(x; t) 



p(x; t) = R (x; t) ; p(x; t) 



dx 



(17) 



(18) 



(19) 



(20) 
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we get 



#( ar t)—\ d8( > X '' t) + 
K[X,t) dx[ dt + 

| 1 f ds( X] t) \ 2 | h 2 d 2 R(x;t) 



dt 

-0, (21) 



2m \ dx J 2mR(x; t) dx 2 

which is, as we have already seen|IJ, together with the equation (|3p, equiva- 
lent to the Schrodinger equation 

% 2 d 2 ib(x;t) Tr . , . ,.dib(x:t) 

with 

^){x- 1) = R(x; t)e is{x > t)/h , (23) 

as we wanted to derive. 

One may note that the second expression in the equation (p0|) is an- 
other restriction, since we are imposing that, in our equilibrium situation, 
the macroscopic momentum has to be represented by the gradient of some 
function. 

From what we have said above we may write the new axioms of the theory 

as: 

Al: Newtonian particle mechanics is valid for every individual system com- 
posing an ensemble; 

A2: The Liouville equation is valid for the description of the ensemble be- 
havior; 

A3: In a thermodynamic equilibrium situation one must have the restriction 
flTFD valid. Besides, the restriction represented by ( |20|) must be also 
applicable at this equilibrium situation. 

In the next section we will show how these postulates are connected with 
the ones used in the previous derivation M. 
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3 Connection with Previous Derivation 



Within the perspective of our previous derivation |T| we define the infinitesi- 
mal Wigner-Moyal transformation 



p eq {x ± 5x/2; t) = Z Q (x, 5x/2; t) = e ip5x/n F{x,p; t)dp 



(24) 



where we called Zq the characteristic function of the probability density 
F(x,p;t). 



Clearly 



p(x;t) = lim Zq(x,Sx/2]t) 

8x^0 



and 



p(x; t)p(x; t) = lim -h 

ox— >0 



dZ Q (x,Sx/2;t) 
d(8x) 



Also 



p 2 F(x, P] t)dp = lim _ h 2d 2 Z Q (x,5x/2;t) 



Sx^o d(5x) 2 
Then, the right hand side of equation (^) may be written as 



(25) 
(26) 

(27) 



J [p - p(x;t)] 2 F(x,p;t)dp 



lim 

5x^0 



h2 d 2 Z Q (x,5x/2;t) | h 2 fdZ Q (x,5x/2;t)\ 2 ' 



d(5x) 2 p(x;t) \ d(5x) J 

which may be rearranged, after some mathematical calculations, as 



(5p) 2 p(x;t) = —% p(x;t) lim 



d 2 



Sx^o d(5x) 



\n[Z Q (x,Sx/2;t)} 



(28) 



(29) 



Performing the expansion of Zq(x, 5x/2; t) in terms of the density, as 
given in (0), and taking the indicated limit in (p9|), we may write this last 
expression, using @ as 



(Sp) 2 p(x;t) = -—p(x;t) 



d 2 ln[p(g;t)] 

dx 2 



(30) 



which is equivalent to the equation (|T9"P. 
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We note that the first two postulates of both derivations are identical M. 
This means that the formal connection between them has to be settled upon 
their third postulates. We thus see that the third postulate we present in the 
last section is formally equivalent to postulate the adequacy of the Infinitesi- 
mal Wigner-Moyal Transformation together with the restriction, made in the 



previous derivation, of writing the characteristic function (23) as the product 



Z Q (x, Sx/2; t) =4j\x- Sx/2; t)iji{x + Sx/2; t). (31) 

In the appendix we show how the connection between both derivations 
works by means of a concrete example. It will be seen that the restriction 
(|3~TD, together with the first equality between p eq and Zq, is equivalent to the 
impositions made in the third postulate presented in the last section. 

It is also interesting to note that we may consider the characteristic func- 
tion Zq(x,5x/2;t) as a momentum partition function, in the same mathe- 
matical sense that the function Z of usual statistical mechanics is an energy 
partition function. 

Indeed, we may establish the close formal use of these two functions for the 
obtainment of useful statistical quantities. Such a comparison is represented 
in Table I, where we present both the cases Z and Zq in contiguous columns. 



4 Statistical Dynamics 

From what we have already said we are now in position to obtain an inter- 
esting equation expressing the dynamical behavior of the system statistics, 



as considered in terms of the functions x and p(x; t) above defined 118 



To achieve this goal we begin considering the mean statistical energy 



H(x; t)p(x; t) 



p 

2m 



+ V(x) 



F(x,p; t)dp. 



(32) 



where the mean is taken over the momentum space only. Using the relations 
©j (H) and ( PI ) we may write 



H(x;t) 



p{x;tf 
2m 



+ V{x) 



fr 2 d 2 In p(x; t) 



■5m 



dx 2 



(33) 
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In this case, the Liouville equation @ gives 

dp(x; t) d 
dt dx 

which may be compared with the Hamilton equation related with the indi- 
vidual systems. We then may call the second term on the right the statistical 
potential and our calculations show that it is related with the momentum (or 
energy) fluctuations. It is precisely this statistical potential that will bring 
the systems to (stable or metastable) thermodynamic equilibrium. 

However, we shall once more stress that his equation reflects the dynamics 
of the statistical variable p(x;t), since the individual ensemble constituents 
satisfy Newton's equations, as postulated. Indeed, for a dispersion-free en- 
semble we have p(x; t) = p. This means that the dispersion (|8|) is identically 
zero and we regain Newton equation (in Hamilton form) from fl34|). We have 
already shown this property in our first paper |]] by other ways. 



H(x;t) 



d 2 p(x; t) 



lmp(x; t) dx 2 



(34) 



5 Metastability: Qualitative 

In the second section we found a relation between entropy and the momentum 
mean square deviation. This relation was based on the equation ( |T0D where 
the usual connection is made between the entropy and the probability density. 

This relation means, because of the process of derivation, that all the 
quantum mechanical solutions we get for a particular problem, the states, are 
related with thermodynamic stable or metastable equilibrium situations — 
whether we are in an absolute or relative entropy maximum. We expect that 
the excited states are related with metastable situations while the ground 
state represents stable thermodynamic equilibrium. 

To see this qualitatively we first write 

S n (x;t) = k\n[p n (x;t)], (35) 

where n labels some specific state. It is important to stress that each partic- 
ular state will have an associated entropy function. 
We may thus define the Gibbs entropy function as 

G n = p n (x;t)ln[p n (x;t)]dx, (36) 
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where here G n will be more properly a constant, since the probability den- 
sities are time independent — as they are related with equilibrium situations. 
This constant, however, will have different values depending on the specific 
state under consideration. We thus expect that its values are such that the 
Boltzmann theorem applies and the more excited the state is the most neg- 
ative the value of G n — for spontaneous transitions between these states will 
imply in the decrease of this function values. This last definition is sufficient 
to take a glimpse on the process of met ast ability. 

As an application we may consider the harmonic oscillator problem. Thus, 
we may use the probability densities related with its states and plot the 
corresponding values of the Gibbs entropy function. The result is given in 
table II. We note that the metastability condition 

AG = G n - G m < 0, (37) 

is satisfied when n > m, n > 0, implying that these states are thermody- 
namic metastable ones. If n = we obviously have a stable thermodynamic 
equilibrium. We also see, from table II , that the more excited a state is the 
smaller the difference AG with respect to its neighbor n±l states. This may 
be interpreted as giving the increase in the probability of transition induced 
by the fluctuations. It is also interesting to note that the values of table II 
have a exponential-decay behavior, as expected. 

This qualitative analysis thus suffices to have an account of metastability. 

6 Metastability: Quantitative 

We may deal with the metastability problem directly from the three fun- 
damental postulates with the only modification that, when metastability 
is being taken into account, the Liouville equation has to be replaced by a 
Master equation. We will use the previous derivation method throughout [|l| 
because it is the most adequate for formal calculations. 

Thus, we begin by considering that the system is in one of its metastable 
equilibrium states, labeled by i, and that it is able to make transitions to 
some final states /. Since we are talking about spontaneous transitions, we 
impose that the energy of the initial state be greater than those of the final 
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ones. The Liouville equation becomes 

= -£*(«'-/)*; (38) 

where 

ei > e f . (39) 

Using the same definition of the characteristic function (|24|) and following 
the same steps of the demonstration previously done0, we arrive at 

at m oxo(ox) ox r 

where we put 

i = £i2(i-/). (41) 

Using again [JXJ] the ansatz 

Z Q (x, 5x/2; t) = ct)\x - Sx/2; t)<f>{x + 5x/2; t), (42) 

writing <f>(x;t) as in (|23"D, expanding (f42|) up to second order in 5x and sub- 
stituting this expansion into equation fl40|) we get, in zeroth and first order 
in &c, the two equations 

dl? + d_ ( R2 ds/dx\ 
dt dx \ m ) t 

and 

dt 2m \ dx J V } 2mR{x; t) dx 2 ' 1 ' 



We note that equation (f43|) represents a continuity equation with a sink. 
This equation can be cast into a more interesting appearance if we put 

R(x;t) = i?i(x;t)e~* /T , (45) 

for, in this case, equation (H3|) may be rewritten as 



ii 



while the equation (pH]) remains the same with the function Ri(x;t) instead 
of R(x; t) 

ds{x-t) J_ ( ds(x;t) \ 2 _ h 2 d 2 R 1 {x-t) 

dt 2m { dx J { 1 2mRi(x;t) dx 2 { ! 

Equations (^7)) and (pRf ) are equivalent to the Schrondiger equation 

h 2 d 2 ip n {x;t) dip n {x;t) 

~ 2^-dx 2 — + V(X)MX] t] = ST - ' (48) 

if we put 

M^t) = Ri(x;t)e ls ^\ (49) 

and represent the (metastable) thermodynamic equilibrium of the system. 
The complete solution to the problem is, thus, 

4> n {x-t) = Mx;t)e~ t/T , (50) 



with the mean lifetime r defined as in fl4ip. This is the usual quantum 
mechanical result found by other means. 



7 Conclusion 

In this paper we have shown how the previous derivation of the Schrodinger 
equation may be accomplished by more sounded physical principles, which 
unraveled some of the characteristics of the quantum formalism — mainly its 
relation with thermodynamic equilibrium calculations. 

The appropriate epistemology underlying all those achievements will be 
postponed to a future paper. We hope, however, that this second derivation 
process makes the foundational principles of the theory most appealing. 



A Boltzmann Equilibrium 

We want in this appendix to apply some of the concepts developed in the 
main text to the special case of a Boltzmann distribution. We will use the 
Infinitesimal Transformation approach]!]], since it is more straightforward, 
when formalism is being considered. 
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We thus consider an ensemble of systems (S) in contact with a reservoir 
(0), called the heat bath, responsible to maintaining the temperature of (S) 
fixed on some value T. The interaction is considered sufficiently feeble to 
make it possible to still attach a hamiltonian function for (S), not depending 
upon the degrees of freedom of (O). In this case, in the thermal equilibrium, 
we have the canonical probability distribution function given as 

F(q,p) = Ce- 2 ^ p \ (51) 

where H(q,p) is the hamiltonian describing the ensemble and 

2/*=^, (52) 

with Kb being the Boltzmann constant, T the absolute temperature and C 
some normalization constant. 

The Hamiltonian may be written as 

N 2 

H(q,p) = Y,7r^ + V( qi ,..,q N ), (53) 

n=l lm n 

where we are assuming that the systems composing the ensemble have N 
degrees of freedom, their constituents have masses m n ,n = 1..JV and are 
acted by a potential V(qi, ..,?jv), not depending upon the velocities nor the 
time. 

Using the Infinitesimal Wigner-Moyal Transformation 

Z Q (q,5q/2) = CJ F(q,p)e l ^ pMh dp, (54) 
the characteristic function becomes, after the integration, 

Z Q (q,5q/2) = C 1 e-^ v ^-^e~^^ {5qn) \ (55) 

where C\ is the constant C modified by the integration process. 
This characteristic function is, clearly, a solution to the equation 

" ft 2 SfZ Q L ^dV 
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obtained|ffl| from the Liouville equation by means of the transformation (|54|). 

As we have already said|]], it must be possible to write the characteristic 
function as the product (see equation (53)) 



Z Q (q, Sq/2) =^(q~ Sq/2; t)i/>(q + Sq/2; t), 



(57) 



to derive the Schrodinger equation related with the possible thermodynamic 
equilibrium distributions. In this case it is easy to see that the probability 
amplitudes ip must be written as 



(58) 



i/;(q-t) = Jc ie -^ v ^'-^e- iEt/h , 



giving, for the characteristic function (up to second order as used in the 
Infinitesimal Transformation derivation^), 



Z Q (q,Sq/2) = C 1 e 



-2/3 



K(gi,..,9Jv)+5E(^«) 



2 a 2 v 

d 9n 



(59) 



Comparing the expression (^) with (|58|), we see that the possibility of 
writing the characteristic function as the product (|57|) is related with hav- 
ing, in the vicinity of the point q = (q 1} ..,q n ), where this function is being 
calculated, 



d 2 V 



<5(Jn=0 



But the expression (|59|) is equivalent to take 



Pe q {q) 



-2pV(qi,..,q n ) 



(60) 



(61) 



as the probability density in configuration space and write the characteristic 
function as 

Z Q {q, Sq/2) = Peq {q± Sq/2), (62) 
representing the first equality in (El), if, and only if, 



dV 



dq n 



0. 



(63) 



Sq„=0 



in the considered point. 

The point where conditions ( |6~0"1) and (|63|) apply defines a mechanical equi- 
librium point for the considered systems. That is, the characteristic function 
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for this problem, where we consider an ensemble of systems (S) in ther- 
mal equilibrium with a reservoir (0), is equivalent to the probability density 
function taken on points slightly apart from the mechanical equilibrium ones. 
The expression ( |62D is just the one we used in section 3, equation (0), in 
the second derivation method. 

Here the connection between the derivations becomes very clear. In our 
specific example the entropy is just 

S( qi ,.,a- N ) = - V{qi - qN \ (64) 

and, since the temperature is fixed, the thermodynamic equilibrium is at- 
tained when ( |5"B"[ ) and ( EDI ) are valid, which means that the entropy is a 
maximum — which is also the mechanic equilibrium situation for this prob- 
lem. We note that the equation ( |5T)D is also equivalent to equation (|H|). 

We may go one step further and ask about the Schrodinger equation re- 
lated with the solution given by the amplitude fl5lf ). In this case, substituting 
this amplitude into the Schrodinger equation 

N h 2 d 2 ti 

- E ^ZTtJ + QnW = E ^ ( 65 ) 



we get 



^ 2m n dql 
*{3h 2 d 2 V MV( "(3 2 h 2 



dq n 



E, (66) 



n=l 2m n d( il f^l 2m n 

which, using the relations ( |60|) and (|63|) , becomes, in the mechanical equilib- 
rium point q = (ql, .., q%), 

E = V(q° 1 ,..,q° N ) + NK B T, (67) 

as expected for this problem. 

The last term on the right is related with the energy given by the reser- 
voir (O) to the system, for each degree of freedom — and is such as if this 
reservoir were constituted of iV identical independent harmonic oscillators 
contributing, by equipartion of the energy, with the amount KbT. If the 
temperature is zero, then it is the same as if the reservoir doesn't exist. In 
this case the energy becomes 

E = V(q° 1 ,..,q° N ), (68) 
as should be in a mechanical equilibrium situation. 
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A.l Harmonic Oscillator 



We may particularize our example and consider the case of an ensemble of 
systems comprising N independent harmonic oscillators. We thus have the 
potential function of these systems given by 

1 N 

V(q u .., q N ) = - m n ulql. (69) 

A n=l 

Equation (|63|) gives the equilibrium points as qi = 0, i = 1,..N, while 
equation fl6T| ) implies that 

K B T = lhu n , (70) 

which means that all the frequencies have to be equal. 
This last equality gives, for the energy, 




(71) 

and, for the momentum mean square deviation, see equation ([18]), 



^ = hu. (72) 

We may, thus, interpret these results as following: the reservoir (O), 
when fixing the systems temperature T, gives to each oscillator the amount 
of energy represented by fl7D|). These oscillators, even in their mechanic 
equilibrium situation will have their momenta and displacements fluctuating, 
because of the energy given by the reservoir (O) to attain thermodynamic 
equilibrium. 

When the temperature approaches absolute zero, the frequencies also 
approaches zero — because of (|T0D — signifying that all oscillations tend to 
stop. In this case the momentum fluctuations also tend to disappear. Re- 
turning to equation ( |i4|) we see that there will be a density fluctuation as 



small as the temperature is. The probability of having a mean quadratic 
displacement obviously also tends to zero. In this case, having the fluctua- 
tions in the momenta and the displacements both tending to zero violates 
the thermodynamic equilibrium condition flT7D, and we conclude that, at 
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this temperature, the oscillators are in mechanic equilibrium but are not in 
thermodynamic equilibrium — if they were to be represented by a Boltzmann 
distribution function at such a temperature, which is not the case. Indeed, 
this last odd result may serve as an indication of the failure of the ensemble 
statistical description by means of the Boltzmann distribution at such low 
temperatures. 
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Energy Partition Function 


Momentum "Partition Function" 




Zq = f e ipdx / n F(x,p;t)dp 


^ ~ df3 


p - hm Sx ^ ih d( ( 5x f 


~P2 1 d 2 Z 

^ ~ Z d/3 2 


_ i lTT1 h 2 9<z Q 


(A£) 2 - 9 JP 


(5 P r = \im 5x ^-h 2d -^ 



Table 1: Comparison between the energy partition function as usually defined 
and the momentum "partition" function as defined in the present work. 



level (n) 


Gibbs Entropy (G n ) 





-1.07237 


1 


-1.34273 


2 


-1.49859 


3 


-1.60978 


4 


-1.69650 


5 


-1.76803 


6 


-1.82901 


7 


-1.88216 


8 


-1.92927 


9 


-1.97179 


10 


-2.01020 



Table 2: Gibbs Entropy values for the harmonic oscillator. 
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